Quantum dynamics, dissipation, and asymmetry effects in quantum dot arrays 
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We study the role of dissipation and structural defects on the time evolution of quantum dot 
arrays with mobile charges under external driving fields. These structures, proposed as quantum 
dot cellular automata, exhibit interesting quantum dynamics which we describe in terms of equations 
of motion for the density matrix. Using an open system approach, we study the role of asymmetries 
and the microscopic electron-phonon interaction on the general dynamical behavior of the charge 
distribution (polarization) of such systems. We find that the system response to the driving field is 
improved at low temperatures (and/or weak phonon coupling), before deteriorating as temperature 
and asymmetry increase. In addition to the study of the time evolution of polarization, we explore 
the linear entropy of the system in order to gain further insights into the competition between 
coherent evolution and dissipative processes. 

PACS numbers: 73.21. La, 73.63.Kv, 85.35.Gv 
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I. INTRODUCTION 

The size reduction of devices in microelectronics, and 
the fundamental limitations encountered or anticipated 
due to quantum mechanical effects, have turned attention 
to understanding nanometer scale structures. Systems 
where intrinsic quantum mechanical effects can play an 
important role and could even be used to advantage are 
under intense scrutiny. One recent prominent example of 
such systems are quantum dot cellular automata. 

The original concept of "cellular automata" was in- 
troduced as n-dimensional arrays containing finite am- 
plitude (information) per cell (or site) and connected to 
one another according to certain rules. This simple defi- 
nition, however, governs-|the evolution of the entire array 
in rather complex ways.EJ Quantum cellular automata are 
a natural generalization of this concept, where each site, 
contains a quantum mecharujcal probability amplitude.EI 
In 1993, Lent and coworkers,cl proposed the use of cellular 
automata architectures, composed of nanometer-scaled 
quantum devices (quantum dots) coupled through care- 
fully chosen hopping and Coulomb interactions, to en- 
code classical binary information in the different charge 
configurations of the system. The typical basic element 
in these quantum-dot cellular automata (QCA) is a cell 
consisting of four quantum dots (QD) located at the ver- 
tices of a square and connected via tunneling barriers to 
their neighbors (see Fig. |l|a) . For two mobile electrons in 
each four-dot cell, Coulomb repulsion between the elec- 
trons causes the charge in the cell, in the ground state, 
to align along either of the two diagonals. The symmet- 
rical arrangement of the system means that these two 
"polarization" states are degenerate and can be used to 
represent logic and 1 as the bits in this system. The de- 
generacy is split by an appropriately designed driver field 
which allows control of the polarization of the ground 



state (and typically implemented as a second identical 
cell with externally controlled polarization). Tunneling 
barriers between dots are designed in such a way that 
intra-cell tunneling is possible, but mter-ccU tunneling 
is not, and only Coulomb interaction from cell to cell is 
possible. Elemental cell designs have been experimen- 
tally realized using metal islands, and logic operations 
AND, OR and NOT have been imnLemented using ar- 
rays of such cells suitably arranged.QtrQ Also, a clever im- 
plementation of QCA using magnetic elements has been 
recently reportedE 

Although originally proposed as 'classical' and near 
dissipationless bit operators, QCA could in principle also 
be used tp implement a quantum computer, as discussed 
recentlyJ3 if the entire array system is coherent enough 
(and for sufficiently long times). In either the classical 
or quantum operation mode, QCA computers (as well as 
any other) should keep errors to a manageable limit. Of 
special significance in these systems, decoherence and dis- 
sipation are external sources of 'error' originating in the 
coupling between the system and its surroundings. Other 
aspects that have to be considered as possible sources of 
loss of control or desirable operation are imperfections 
in the dot fabrication, stray charged impurities, and the 
role of switching fields or drivers. It is then of interest 
to study the dynamical evolution of the QCA and mon- 
itor the charge density distribution over the basic cell 
elements under the influence of a driver cell, with empha- 
sis on the effects of imperfections and dissipation. The 
charge configuration iS|,conveniently monitored through 
Lent's cell polarization^ which is essentially a measure of 
the degree of alignment of the charges along the diagonals 
of the square cell: 
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FIG. 1: (a) Energy level structure vs. t for a target cell with 
driver at +1. Energy in units of 1/ ~ 1 meV, and 5 = 0. 
Sketch at bottom: model QCA with four quantum dots per 
cell. Two charges are allowed to tunnel between nearest neigh- 
bors, with probability t. Driver cell on right indicates a state 
with polarization Pdriver = +1. (b) Level structure vs. dot 
size asymmetry 5. Notice anticrossings for (5 ~ 0.1 and 0.2 in 
excited states. Separation to driver, c = a in both panels. 



where pi is the charge density at site i (see Fig. |l^), 
and the numbering of the sites is in clockwise order. 
In previous works, the effect of dissipation on QCA has 
been included in the quantum equations of motion in a 
phenomenological way,u through a damping parameter. 
This model is shown to give qualitatively correct results 
only in the low temperature limit, as it will be evident 
in our calculations. Also, the temperature dependence 
of the cell polarization has been studied for the case pf 
a cell composed of two Coulomb-coupled double dots,t2l 
and the contribution of relaxation processes to the speed 
of the cell response was discussed using stochastic equa- 
tions for the population differences and dipole moments 
of the double dots. Thcj^ffieci^ of imperfectkms have been 
considered before, withtillljllj and withoutliil dissipation, 
and found to have strong effects in the dynamical re- 
sponse, by delaying or slowing down the response of the 
target cell, and producing a less than optimal desired 



polarization. In this paper, we include the important ef- 
fect of the environment through a microscopic model of 
electron-phonon interactions. The thermal bath is mod- 
elled with phonons in equilibrium at a given temperature 
and structural asymmetries are introduced via size vari- 
ations of the quantum dots. 

We should point out that in addition to their relevance 
as a new computational paradigm, QCA arrays provide 
a conceptually simple, yet intricate and subtle, quantum 
mechanical system in which to study fundamental coher- 
ence and dissipative processes. The parametrical require- 
ments to make practical QCA devices yield a finite set 
of two-particle states which can be manipulated to yield 
insights into the loss of coherence in these and general 
quantum systems. As such, QCA provide then a model 
condensed matter system where dynamics and decoher- 
ence can be explored in detail. To-jtliiA-jeffect, we use 
a quantum open systems approachll3il3ila to study the 
dynamical evolution of the reduced density matrix ps{t) 
of the system, and evaluate the competition between co- 
herent (driver cell) and incoherent (thermal bath) field 
effects in the response of the basic cell. We study the ef- 
fects of temperature, electron-phonon coupling intensity, 
and structural asymmetry on the behavior of the polar- 
ization as a function of time. In particular, we discuss 
the asymptotic polarization after switching events, and 
once the driver is in a stationary configuration, for dif- 
ferent temperatures and structure parameters. An anal- 
ysis of asymptotic fixed points in the resulting dynamics 
provides information on the system, as we will explicitly 
demonstrate and discuss later on. 

A measure of the mixed character of a system described 
by a density matrix p is provided by the linear entropy^ 
s[p\ — Tr{p — p^^Ln This quantity, which can be thought 
of as the lowest order approximation to the von Neumann 
entropy, is also called the idempotency defect or the pu- 
rity of a state. Notice that the linear entropy is zero for a 
pure state, since Tr{/5} = Tr{/5^} = 1 in that case. Non- 
zero values of s[p] provide then a quantitative measure of 
non-purity of the state of a system. When monitored in 
time, the linear entropy provides a convenient measure 
of how fast the loss of quantum purity occurs in a system 
in contact with a bath, or the increase in coherence un- 
der appropriate 'pumping' conditions, as we will discuss 
below. Other uses of s[p] include the study of quantum 
manifestations of chaotic behavior in a complex system.tS 
We explore the influence of temperature, electron-phonon 
coupling intensity and driver cell switching times on s[p\, 
in order to understand the process of dissipation in the 
QCA system. 

We find in all cases that high dissipation and asym- 
metry influence adversely the response of the basic cell, 
as one would intuitively anticipate. However, the behav- 
ior with temperature is quite subtle. For example, low 
temperature results in better switching and more faithful 
following of the driver by the target cell. Higher temper- 
atures eventually produce quenching of the polarization 
and a dynamical evolution which basically ignores the 
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driver. On the other hand, size defects in the dots result 
in lagging of the switching behavior with respect to the 
driver, and a faster switching window. This is produced 
by the pinning of a given polarization state by the ener- 
getics of the size defect, so that the driver potential has 
to apply a strong force before the target cell can follow. 
Similar lag and/or pinning, as well as sudden switch of 
the polarization is seen in the case of two target cells be- 
ing driven by a single cell. A given cell in the middle 
of a chain would clearly have to deal with the influence 
of more than one cell, resulting in a more complex po- 
tential landscape and corresponding dynamical behavior. 
We also find for a single cell that in the case of low tem- 
perature or weak dissipation the cell recovers dynamical 
coherence, due to the 'pumping' of the driver, clearly 
overpowering the thermal bath. Other regimes will be 
discussed in detail. 

The behavior of the polarization of the basic cell as 
a function of time when the polarization of the driver 
is fixed, as well as wheruJ.t changes linearly, has been 
briefly discussed recentlyJ13 There, we studied the sym- 
metric (no defect) case, but taking account of dissipation 
in the system. In what follows, we present a full report 
of the QCA dynamical behavior under different electron- 
phonon and temperature conditions, and analyze the lin- 
ear entropy of the system during these processes. 

The paper is organized as follows: in section || we de- 
scribe the model for the entire system, composed of the 
QCA and the thermal bath. We use an extended Hub- 
bard Hamiltonian to describe the Coulomb interactions 
present in the system (target cell(s) -I- driver cell). The 
thermal bath is represented by a set of harmonic oscilla- 
tors, and the electron-phonon interaction enters as emis- 
sion and absorption of phonons during electron tunneling 
events between quantum dots in the system. A Marko- 
vian approach to the dynamic al evo lution of the reduced 
density matrix is presented in |II A| , where the transition 
rates for the electron-phonon interaction in this system 
are explicitly obtained. In section ^ we present and dis- 
cuss the behavior of the polarization and linear entropy, 
as function of relevant parameters, with special emphasis 
on effects of dissipation and asymmetry. Here, we also 
discuss the dynamics of multiple cells, as the QCA ul- 
timate devices involve multiple cell geometries. Finally, 
we summarize and discuss our results in section [V. 



terms, as well as intra-cell tunneling, 



II. MODEL 

Consider a problem described by the Hamiltonian H — 
Hs + Hji + VsR, where the subscripts S and R stand for 
system and reservoir, respectively. In our case, Hs de- 
scribes the target cell in the presence of a driver cell which 
provides external driving fields on the system. The inter- 
actions are introduced via an extended Hubbard Hamil- 
tonian with both intra- and inter-cell Coulomb repulsion 



Hs = Y^ejfij - J2 {tijc\^j + c.c.) 

j <ij> 



V 1j I L2,l Lj 



(2) 



The first three terms describe the target or active cell and 
the last term describes the Coulomb interaction between 
the electrons in the target cell and a driver cell which 
is externally controlled. Here, is the confinement en- 
ergy level in the j*'* quantum dot (QD) in the target cell. 
In this term, possible dot asymmetries or local poten- 
tial variations are introduced by changing the energy in 
a given k*^ QD, — > e — (5, where 6 is the measure of the 
imperfection due to different dot size or a change in the 
local environment, and e is the energy of the other (iden- 
tical) dots. The amplitude is the tunneling matrix 
element between nearest neighbors < ij > on the same 
cell, c| and cj are creation and annihilation operators for 

site j, hj = CjCj is the electron number operator at site 
j, Vij is the Coulomb interaction between sites i and j 
in the target cell, and Wij is that between site i in the 
driver cell (with charge density nf{t), which in general 
changes in time) and site j in the target cell. A gen- 
eralization to multiple-cells is straightforwardly given by 
inter-cell interactions.L3 Following the notation of Lent et 
a/.,B the distance between nearest-neighbor dots (length 
of the square side) is a, while the separation between the 
target and driver cells, located side by side, is c (see inset 
Fig. la). Notice we consider spinless electrons, although 
considering two-particle states with given total spin is 
also possible in our formalism. 

The Hamiltonian for the reservoir is given by a set of 
harmonic oscillators of frequency ujk, Hr — f>^kb\bk^ 
where S|,(6fe) are creation (annihilation) operators of 
phonons. Finally, the electron-phonon interaction is 
given by a general expression. 



SR 



2^ 2^ 0^k^jC\Cj[bl 
k <ij> 



(3) 



corresponding to emission/absorption of a phonon for 
electron tunneling events and with matrix elements given 
by the coefficients akij- Notice that the physics of the 
phonon interactions in the system is well described by 
the model chosen here. Although phonons are in prin- 
ciple present everywhere, and not only during tunneling 
events, they only provide an 'on-site' energy renormal- 
ization/shift and broadening, which will be ignored here, 
as it yields a uniform correction to the entire cell. In 
contrast, the phonon- assisted tunneling events described 
by VsB. affect significantly the dynamics of the cell, as 
they provide an effective coupling mechanism among the 
different levels in the spectrum and allow for energy re- 
laxation processes. 

We assume, as in ||l^, that the quantum dots have a 
characteristic size d ~ 50 nm defined on a GaAs / AlGaAs 
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heterostructure, with effective mass m* = 0.067. The 
typical distance between nearest-neighbor dots is a > 
100 nm, and the dielectric constant of the medium is w 
12. These values give an estimate of the Coulomb repul- 
sion between nearest neighbor dots oiV \ meV, which 
is taken as the unit of energy in this paper. One should 
notice that the first excited single-particle state in each 
dot hes at quite a high energy, A « (?i^/2m*)(27r/d)^ w 
9 meV ^V, so that one can safely consider a single or- 
bital per site, although additional single-particle levels in 
each dot could in principle be included. Notice also that 
the on-site Coulomb energy is quite high and prevents 
the double occupation of the dot, so that our basis set 
of states ignores that possibility. We should comment 
that successful implementations of QCA elements using 
metallic islands are in a multiple-orbital regime, since in 
that case IS. <^V for typical l/im-size islands. The con- 
dition of overall electrical neutrality of the target (and 
driver) cell is included in the calculation by distributing 
a uniform positive background charge -\-2e throughout 
the cell (e/2 in each dot), which results in multipolar 
fields from/at each cell. 

Notice that we study the behavior of the target cell 
interacting with the driver via Coulomb interaction 
only (no inter-cell tunneling allowed). This interaction 
is described by the last term in the Hamiltonian 
Eq. (^), and assumed to be time dependent via nf{t). 
For a fixed driver and no dissipation, the solution can 
be found by direct numerical diagonalization of Hs- 
In that case, the effect of various parameters on the 
polarization have been studied to show that the energy 
spectrum is strongly influenced by even small defects, 
which affect the bistability of the cell.Ey In the case of 
a time-dependent driver, we have to take into account 
that the driver promotes transitions to intermediate 
excited states. In fact, the dynamical competition 
between the driver (coherent terms in the equation 
of motion) and the bath (dissipative terms) causes a 
non-trivial time evolution of the system. We have stud- 
ied different 'driving schemes' but will mostly consider 
here the case of a driver with a linear switching of its 
polarization, Pdriver{i) = 1 ^ 2t/T, for times < t < r, 
where r is the switching time. This corresponds to the 
configuration of charge densities in the driver cell of 
ni{t) = ni{t) = l~ t/r and ni{t) = nfit) = t/r. For 



for nearest-neighbor dots, while all others scale with 
their (inverse) separation. We also set all tij equal to 
a constant t between nearest neighbors only (not to be 
confused with time). 

A. Markovian master equation 

In the interaction picture (/), the evolution of the total 
density matrix is given by 



dt 



1 



4 / dt' 



where ^he cell-bath interaction operator, and the 

expansion is valid up to second order in Vj^^ Jll 

Our approach uses a closed and time-local equa- 
tion for the reduced density matrix (RDM): Pg\t) = 
Trfl{/5'^^^(i)} where the trace is carried out over the de- 
grees of freedom of the reservoir, R. The following fun- 
damental assumptions are made in this approach: (1) 
the system S and reservoir R are initially uncorrelated, 
which imphes p(0) = ps(0)/5j^(0) = /5^^^(0); (2) a station- 
ary reservoir at temperature T with equilibrium density 
matrix pr{Q) = exp(-iJfl,/fcBT)/Tr{exp(-iJii/fcsT)}, 
exists at all times; and (3) the correlation time in 
the reservoir is much shorter than characteristic times 
for the RDM system to change appreciably, meaning that 
p"g\t') K, ffp {t), if t' —t < Tc, while the state of the reser- 
voir at time t' has already no correlation with the state 
at time t (Markov approximation). Thus one arrives at 
the time evolution equation for the RDM element^Cd£- 
tails of the derivation can be found in the literaturetSEj) , 
which in the Schrodinger picture is given by 



'Ps{t)s 



'Ps{t)s 



Rs 



•.Ps{i)r. 



(5) 



the four-dot cell arranged on a square, we take Vij 



V 



The first term on the right represents reversible dynamics 
("coherent" effects), in terms of the transition frequen- 
cies uJss' — {Eg — Eg/ ) / h, while the second term describes 
relaxation (irreversible dynamics). In this part, the re- 
laxation tensor Rss'mn can be written explicitly as 



y -Iss' OmsOns', (s 7^ S ) 



Here, Wnm are transition rates from state \m) to can be expressed as 

which in terms of the properties of bath and target cell, 2,]- 



Wmn = ■^£'^|.g(t^mn)P^^(Wmn) {|'S'„.mpn(a;„„) 

\S,^n\^\\+n(lO„,n)\) (7) 
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where Smk are matrix elements of the electronic part 
of the target cell interacting with the phonon reservoir, 
Eq. (|); n{uj) = (St(u;)S(w)) = (e'""'" - l)-^ is the av- 
erage phonon number in the reservoir at temperature 
T = l/ksfi. The first term (ex n) describes the absorp- 
tion of phonons, while the second (oc n + 1) describes 
phonon emission, both at a frequency given by the sys- 
tem's transition energy ujmn- Notice that the matrix ele- 
ments Smn embody the appropriate selection rules for the 
different states involved in the transitions. For simplicity, 
we consider a model where the amplitude of interaction 
does not depend on the electronic states of the target 
cell, and has the form akij = Dg^^ where D is the con- 
stant prefactor of the electron-phonon interaction jEiid is 
determined from the deformation potential modelH We 
should emphasize that this model introduces dissipation 
effects only at finite temperatpje, and cannot account 
for elastic decoherence effects.l23 Here, gu = g{uJk) ~ 

1/2 

ujfj describes the frequency dependent amplitude of the 
electron-phonon interaction in the deformation potential 
model; T>{ijj) ~ is the density of states in a Debye 
model and oJmn is the transition frequency between states 
m) and |n). It is clear that other phonon channels with 
different dispersion relations and/or matrix elements ex- 
ist, but we believe that other modes would not be as 
effective in introducing dissipation in this system (e.g., 
optical phonons occur at energies too high to couple ef- 
fectively), nor qualitatively change the results discussed 
here. In Eq. (^), 7^5/ is a so-called non-adiabatic parame- 
ter whose real part gives a contribution to the time decay 
of the off-diagonal density matrix elements, and is then 
directly responsible for the loss of coherence. These pa- 
rameters can be written in terms of the transition rates by 
Re -iss' = {Y^k^s ^fcs + Y^k^s' Wks')l'2.- The imaginary 
part of 7ss' is an intrinsic relaxation rate for each transi- 
tion, which we assume negligible, Im ^ss' = 0, as itpeould 
be incorporated in the transition frequencies, Wss'.tj 

It is interesting to point out that Eq. (|^), written for 
the diagonal elements, takes the form: 

'Ps{t)ss = E Ps{t)mra " Ps{t)ss E ' 

m m 

which has a simple interpretation: the probability of find- 
ing the level \s) occupied at time t, ps{t)ss^ increases due 
to transitions from all other levels |m) to \s) (first term) 
and decreases due to transitions from \s) to all other lev- 
els |m) (second term). 

Note also that the condition Ps{t)ss' = in Eq. (|^) cor- 
responds to 'fixed points' of the dynamics of the system, 
as both the terms related to coherence, {ps)ss' for s ^ s', 
as well as the diagonal terms which determine the popu- 
lation in each state, do not change in time. These fixed 
point values depend only on the properties of the transi- 
tion rates and the energy spectrum of the system (as the 
external drive reaches a stationary configuration), and 
provide a direct determination of the dissipation times 
produced by the thermal bath. It can be seen from Eq. 



that the transition rates satisfy 

nra \ /^\ 
Wmn \ ksT J 

which guarantees the detailed balance condition. It is 
important to emphasize that we have derived an explicit 
expression for the transition rates Wmn , which would be- 
have at low temperatures, T — > 0, and/or low dissipation 
(small D^) as the form used by Zozoulenko et ai, in a 
qualitative way.EJ Notice however that the normalization 
condition in pT[ is an added requirement in their phe- 
nomenological approach, while it is included naturally in 
our treatment. We should also mention that the detailed 
balance ratio in (^ will determine the fixed point values 
of ps and all the corresponding physical properties, in- 
cluding the different state populations. Since these will 
in general be asymmetric, the polarization value of the 
QCA is reduced, and as we will see later, this results 
also in a limiting value for the linear entropy. A detailed 
analysis of the asymptotic fixed points in terms of the 
structure parameters will also shed light into the dynam- 
ical evolution of the system, as we will discuss. 

The solution of the master equation (||) for the time 
evolution of the RDM is carried out numerically using 
a Runge-Kutta fourth order algorithm. We ensure, by 
a proper choice of the integration step, that the general 
normalization property Trjps} = 1 for the RDM is satis- 
fied at all times. Notice that an instantaneous basis is in- 
corporated in the calculation to take into account the evo- 
lution of the ground state induced by the time-dependent 
driver configuration. The basis functions are then defined 
by the eigenvalue/eigenvector problem of the coherent 
Hamiltonian of the target cell, Hs{t)\m) = E^{t)\m), for 
a given time t while the driver is in the process of switch- 
ing, and the instantaneous energy in the target cell is 
E^{t). Notice also that the characteristic time variation 
of the driver r is long compared with the target-cell's 
natural frequencies loui, so that t 3> h/ujai. This is effec- 
tively the requirement for an adiabatic switching of the 
driver. As is clear in Fig. |l], the characteristic frequencies 
of the system are uuw ^ 1 meV, so that h/ojui > 1 ps. 
This would require r to be a few nanoseconds for the 
adiabatic switching regime to be valid. 

The numerical solution of the RDM is used to calcu- 
late the dynamics of the polarization in the target cell, as 
the physical observable of interest here. The dynamical 
behavior of P{t), calculated from the charge density in 
each site as pi{t) = Tr{ps{t)ni}, and of the linear entropy 
s[(05], are then determined for a given set of parameters. 
In order to get better insights in the problem, we also 
carry out a calculation of the fully coherent dynamics, 
ignoring the phonon thermal bath (setting D = in Eq. 
[7|). In the coherent limit, the relevant parameters for 
the dynamical evolution are those defining the structure, 
including the tunneling amplitude t, the separation of 
target and driver cell c, the switching time of the linear 
driver r, and the size of the dot asymmetry 6, if any. For 
the dissipative problem, on the other hand, it is clear that 
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in addition to the structure features, the important pa- 
rameters are the temperature of the thermal bath T , and 
the amphtude of the electron-phonon interaction D. Of 
course, the phonon density of states, as well as the inter- 
action details, will also play a role in the evolution of the 
system. Our model parameters provide an adequate and 
realistic description of these systems in semiconducting 
QCA.EI 



III. RESULTS AND DISCUSSION 

Figure l^a shows the level spectrum in the target cell for 
a fixed driver polarization, Pdriver = 1, as shown in the 
inset. We see that the four- fold symmetry of the struc- 
ture and consequent degeneracy is broken by the driver, 
even for i = 0.E3 Moreover, increasing the intra-cell tun- 
neling t produces additional level splittings, as one would 
expect. We should point out that the polarization of the 
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ground state follows that of the driver, 
increasing hopping t reduces that value somewhat. The 
first excited state has J^^T/? — ^^'^ is separated from 
the ground state by a gap of ~ 0.2 in the range shown 
(our energy units are V ^ 1 meV) . The four higher states 
have vanishing polarization, as the charges are located on 
neighboring sites along the sides of the square. The basic 
operation of the QCA requires adiabatic evolution from 
a polarized ground state to the other with opposite po- 
larization. Transitions to higher excited states allow for 
non-adiabatic terms which yield to depolarization of the 
target cell. Figure |l|b shows the influence of a "defect" in 
dot size. Here, the dot at site 1 has an energy lower than 
the others by an amount 5 (corresponding to a larger 
dot). Increasing this asymmetry produces a strong shift 
of the spectrum, and perhaps most important for QCA 
operation, a change of the polarization of the first excited 
state. This is apparent in the anticrossing seen in |^b for 
5 ~ 0.1, where an opposite polarization state takes over 
as the first excited state.E3 The defect then will produce 
a pinning of the polarization when operating the QCA 
if the S is small, and even prevent the switching of the 
target cell altogether if 6 is large. 

When the polarization of the driver cell changes lin- 
early with time, the polarization of the target cell evolves 
as well, as presented in Fig. ^ for the non-dissipative 
case, and for different temperatures. One can see that 
the polarization of the target cell follows the driver very 
well, and switches smoothly over a short time window 
(~ t/5) centered about the time when Pdriver = 0. This 
general behavior is not affected for different values of the 
switching time r we tested (over a wide, range within the 
adiabatic regime; results not shown) .Ej We see that the 
response of the target cell deteriorates for increasing tem- 
perature. At low temperature, T = Ti ~ 1 K, the target 
cell has a slightly decreased polarization after the driver 
is fully switched, and the switching time of the cell is not 
greatly affected. However, for T = 5Ti, the target cell po- 
larization is quenched after 0.5 ns and ignores even a fully 
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lower energy dot, 5 — 0.1, at site 1. In all plots, t = 0.03 and 
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polarized driver. [Notice that the natural unit of temper- 
ature in this problem is the difference in energy between 
the first and second excited states, when the driver is 
fixed at Pdriver = 1 (Fig- 03), as this transition has one 
of the smallest natural frequencies uju' , and it produces 
depolarization of the cell, as discussed above. In partic- 
ular, we take as characteristic its value for t ~ 0.03 in 
the symmetric case, Ti = 0.0862 « 1 K. Temperatures 
are compared to this natural scale in the problem.] The 
deterioration of the response of the cell with temperature 
is consistent with Eq. (|9|) where the different rates at low 
temperature would favor transitions to the ground state, 
while as temperature increases, transitions up and down 
become equally probable. This results in more symmetric 
level populations, more equal electron probabilities in all 
dots, and consequently the depolarization of the target 
cell. It is somewhat surprising that even this low tem- 
perature would yield such strong quenching of the QCA 
polarizability (a value which of course depends also on 
how strong D is in Eq. Figure ^ shows an exam- 
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target cell lagging the switch of the driver increases with S. 
(c) Long-time asymptotic polarization. Pa vs. temperature T 
for different 5 values. In all panels, t — 0.03 and D'^ — 0.05. 



pic of the effects of asymmetry (quantum dot at site 1 
has lower energy than the rest), for the same parameter 
values and temperatures as in Fig. ||a. We observe that 
at low temperatures, the response of the target cell lags 
behind the driver, although the switching time is not 
affected substantially. However, the target cell reaches 
asymptotic values of polarization which are smaller (in 
absolute value) than in the ideal symmetric case. In this 
case too, changing r over a wide range does not affect 
this lagging behavior. It is evident that for T = STi, the 
target cell is almost completely depolarized, despite the 
"help" of the larger dot which pins the polarization of 
the ground and excited states to ~ 4-1 (and induces the 
lag seen at low temperatures). 

Figure ^ shows more of the effects of asymmetry and 
temperature on the response of the target cell. Figure 
^ shows the results for no dissipation, and panel (b) 
those for T = Ti. From these figures, we see again that 
the important effect of asymmetry is to produce a lag or 
delay in the response of the target cell as asymmetry in- 
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creases. On the other hand, the effect of temperature is 
to decrease the asymptotic polarization in absolute value, 
associated with the delocalization of electrons within the 
cell. Notice also in Fig. ^ that for finite T, the polar- 
ization of the target cell continues evolving in time (de- 
teriorating), even after the driver is fully polarized. In 
all cases, the polarization reaches asymptotic values Pa 
(long after the driver has stopped changing) which are 
a manifestation of 'fixed points' of the system dynam- 
ics for the quantum dot populations. These long-time 
values of the polarization depend on temperature (and 
indirectly on structure parameters and energy levels) as 
shown in Fig. |^c, for different asymmetries. The nearly 
full depolarization of the target cell at high temperature 
is evident and anticipated. What is somewhat counterin- 
tuitive is the fact that the depolarization increases with 
asymmetry, as the pinning that yields the lag assists now 
in the depolarization by favoring the population of the 
states with the larger dot. 

Figure ^ shows the effects of tunneling coupling in- 
tensity t on the time evolution of the polarization for 
a symmetric cell system {S — 0). We see that increas- 
ing hopping between dots first introduces a smoother 
transition following the linear driver (r = 2.5 ns) and 
a better polarization. Increasing t however, delocalizes 
the charges in the quantum dots and reduces the polar- 
ization of the target cell, as one can see from the long 
time limits in the figure. In order to look more closely 
at the non-monotonic behavior of Pa, we analyze this 
quantity as a function of the tunneling coupling t in Fig. 

We find that, for T — Ti, Pa attains a minimum 
value at t « 0.02. From this figure we can see that the 
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drop in polarization at small values of t is due to dis- 
sipation processes: We expect that at low t values the 
electrons would tend to be well localized at the corners 
of the cell (ground state), giving P ~ ±1, depending 
on the driver. This is indeed what is observed at low 
temperatures. However, as temperatures increases, de- 
polarization transitions from the ground state to excited 
states become more likely. As t increases substantially, 
the polarization decreases (in absolute value) due to the 
delocalizing effects of both tunneling and dissipation, as 
expected. However, for t < 0.02, the possibility of tran- 
sitions which in fact populate the ground state increases, 
and therefore improve slightly the asymptotic polariza- 
tion of the system. Qualitatively similar results are ob- 
tained in the asymmetric case (Fig. ^), although with 
lower absolute values of the asymptotic polarization in 
this discussed above (Fig. g and g), and a shal- 

lower oscillation. 

We have also studied the dynamics of two target cells 
and driver in line (Fig. There are a number of in- 
teresting features in this case. We notice first that the 
switching for both target cells lags behind the driver con- 



siderably, as the cells switch only after Pdriver ^ 75%. 
Moreover, the switch occurs over a much shorter time 
scale than for the single target cell (~ t/25). This change 
is understandable if one thinks that the target cell clos- 
est to the driver (cell 1) is in essence between two cells 
with different polarization and has its wrong polariza- 
tion aiamewhat pinned (this is behavior also seen in other 
workllj). Increasing temperature to T = Ti makes the 
switching a bit smoother although still with a sudden 
change, and with not as good final polarization. Notice 
in particular that the asymptotic polarization of target 
cell 2 is poorer than that of cell 1. As the tempera- 
ture increases to T = 2Ti, the target cells barely fol- 
low the drive, and the polarization of cell 2 is basically 
zero at long times. It is clear that control of the tun- 
neling process as the "signal" is transmitted down the 
line would be essential in proper operation of QCA ar- 
rays. Such "clocked switching" has been discussed by 
Lent and cojwprkers as a solution to steer cells into bpLter 
complianccEJ Coherent simulations oLsuch steeringJlJ as 
well as experiments in metallic QCAEj verify that this is 
a good operating scheme. 

Knowing the time evolution of the density matrix we 
can also calculate the linear entropy of the QCA sys- 
tem s[ps] , to monitor the degree of non-purity introduced 
during the switching process by the thermal bath. In Fig. 
we show the results for different temperatures, for the 
case of weak and strong electron-phonon interactions, as 
measured by the value of in Eq. (0). As in previous 
cases, the polarization of the driver cell is switched lin- 
early from -1-1 to —1 in a time r — 2.5 ns. In Fig. |^a, 
the case of 'weak' electron-phonon coupling, — 0.05, 
we see that at low temperatures s[ps] increases mono- 
tonically until it reaches a final stable value. It is in- 
teresting to see that as Pdriver = or Pdriver = " 1 (at 
1.25 and 2.5 ns, respectively), one can see kinks in the 
entropy curve, indicating the switch off of driving forces. 
The long time asymptotic value increases (coherence de- 
creases) with increasing temperature, as one would ex- 
pect. However, at still higher temperatures and/or for 
strong electron-phonon coupling (Fig. ^), s[ps] first in- 
creases and then decreases. In these cases, we see a max- 
imum in the linear entropy curve at ~ r/2, when the 
polarization of the driver cell goes through zero. Notice 
in fact that whenever Pdriver — 0, the level spectrum in 
the target cell shows a great deal of degeneracies. As 
the driver turns negative (for time after t/2) the split- 
tings are reinstated and the effect of the driver seems to 
be one of a coherent 'pump' which produces a drop in 
the purity, even as dissipative transitions are present. In 
other words, the thermal dissipation seems to decrease 
with respect to the coherent transitions induced by the 
driver, producing a drop in the linear entropy of the sys- 
tem. After the driver is stationary, after time r, the 
thermal transitions quickly produce the steady state in 
the density matrix (or level populations) given by the 
detailed balance equation (|^). For stronger phonon cou- 
pling, = 0.5, and high temperatures (Fig. |^b), the 
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FIG. 7: Time evolution of linear entropy or purity of the sys- 
tem in a linear switch, r = 2.5 ns at different temperatures, 
(a) Weak coupling, — 0.05, t — 0.03; (b) strong coupling, 
= 0.5, t = 0.03. Notice entropy drop after r/2 at interme- 
diate dissipation. 



purity of the system deteriorates further and the driver 
is not able to 'push' the system back into a more coherent 
state. Notice that the state of 'maximal delocalization' 
would be achieved as psu — ^/n, where n is the size of 
the basis, so that Tr{p^} ~ 1/n. In that case, one would 
obtain s ~ 1 — 1/n, giving 5/6 in our case. This value is 
still higher than the asymptotic value of s for the highest 
temperatures in Fig. so that even then the system is 
not 'fully' decoherent. 

Let us explore this competition between dissipation 
and coherent evolution in more detail. Since Tr{ps} — 1 
is a general property of the density matrix in any rep- 
resentation, one can instead investigate the behavior of 
Tr{/5|} = 1 — sLoj] directly as a function of the rele- 
vant paramcters.Ej Notice in fact that this trace has two 
contributions 

i 

= Y.iPS^^? +Y.\PS^^'' ^ (10) 



namely, the diagonal and non-diagonal elements of the 
density matrix. Comparing the two terms as a function 
of time, for example, would give one further insights of 
the overall complex processes, since the behavior of the 
non-diagonal elements gives information on the evolution 
of coherence in the system. Figure ^ shows the trace as 
function of time for two values of the temperature, explic- 
itly showing the contributions from the diagonal and non- 
diagonal elements of ps- We see in Fig. ^ an increase in 
the contribution from the non-diagonal ( "coherent" ) ele- 
ments, with maximum at t/2, exactly compensated by a 
decrease in the contribution from the diagonal elements, 
as the driver and the bath are promoting transitions. No- 
tice also that the values of the diagonal components be- 
fore and after r/2 are approximately the same (with a 
similar behavior for the non-diagonal terms). Thus for 
low temperatures T — Ti/IQ (panel a), the sum remains 
approximately constant, equal to 1, and the system is 
close to the pure limit with s 0. However, for T = Ti, 
we see that this sum decreases with time, due essentially 
to the corresponding decrease in the diagonal elements. 
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Notice that as the driver switches and 'pumps' the target 
cell system, the off-diagonal coherent terms contribute 
less to the sum as the temperature increases (its peak 
value is smaller in (b) than in (a); so that the system is 
less coherent indeed), and their contribution for longer 
times is \owei than at the start of the switch (time zero). 
This gradual decay in the coherent contribution shows 
even more substantially (and less intuitively) in the di- 
agonal terms, which one would expect to be less affected 
perhaps. Since Tr{/5} = 1, the drop in the diagonal terms 
of Tr{/5^} means that the population of each level {psu) 
has been re-distributed and made somewhat more equal 
(although still much higher than the 'maximal delocal- 
ization' state discussed above with Tr{/5^} ~ 1/6). 

It is also interesting to investigate the influence of 
switching time and of the function controlling the switch- 
ing. Figure ^ shows results for Trj/Sg} for a smoother 
driver transition. In this case, the driver has no rapid 
off/onsets and the evolution of Pdriver is more adiabatic 
near the ends of the sweep than before (see panel c). 
Panel (b) shows the trace at low temperatures, which 
shows a similar rise (drop) in non-diagonal (diagonal) 
density matrix elements as in Fig. However, as the 
time during which Pdriver — is shorter now (effectively 
larger slope in panel c), the window of variation is also re- 
duced. Higher temperature in panel ^ shows also a nar- 
rower time window variation, but an identical amplitude 
as in Fig. ||b. Other runs with different driver switching 
slopes follow this behavior too: narrower windows for 
larger slopes near Pdriver — (not shown). This would 
suggest that the amplitude variation of the diagonal/off- 
diagonal components is only a function of temperature, 
while the time window over which they vary is a function 
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of how long the driver is close to zero polarization. 



IV. CONCLUSIONS 

We have studied the behavior of the polarization of a 
square array of quantum dots with two mobile electrons. 
The quantum dynamical evolution in response to a driver 
cell whose polarization is changing with time in a pre- 
scribed manner, has been shown to depend strongly on 
switching or driver characteristics, symmetry of the cell 
and on temperature. Our results show that the asymp- 
totic values of the polarization in a stationary driver are 
fixed points of the system that change with tempera- 
ture, tunneling probabilities and imperfections, but not 
on driving schemes. The optimal response of the QCA 
arrays is affected adversely by temperature and imper- 
fections, with a strong tendency toward depolarization 
as temperature or asymmetry increase. Our calculations 
indicate that good control of the QCA is only achieved for 
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temperatures of at most a few kelvins. From our two tar- 
get cell results, it is likely that "clocked switching" would 
be required for faithful operation of a chain, so that the 
polarization information is not degraded down the line. 
We have also analyzed dissipation effects on the response 
of the basic cell through calculation of the linear entropy 
or purity of the system. We find that the study of this 
quantity is a good tool in furthering the understanding 
of a quantum system as relatively simple, but still quite 
subtle, as the QCA architecture. Apart from providing 
insights as a function of time and driver characteristics, 
the entropy also shows clearly that at high temperatures, 
the basic cell is no longer in the required/desired final 



state but in one where the polarization tends to zero and 
maximizes the entropy. 
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